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ABSTRACT 

Context. TW Draconis is one of the best studied Algol-type eclipsing binaries. There is significant evidence for miscellaneous physical 
processes between interacting binary components manifesting themselves by period and light curve changes. 

Aims. Obtaining new set of photometric and spectroscopic observations, we analysed them together with the older spectroscopic and 
photometric data to build model of this eclipsing system with respect to observed changes of 0-C diagram and light curve. 
Methods. Reduction of new spectra was carried out in the IRAF and SPEFO programs. Radial velocities were determined manually 
using SPEFO, by CCF and from the program KOREL. Orbital elements were derived with the FOTEL program and via disentangling 
with KOREL. The final combined solution was obtained with the programs PHOEBE and FOTEL. 

Results. Photometry shows small irregularities in light curves as a results of pulsating of one component and spot activity. Using net 
of KOREL outputs we found the mass ratio q=0.405(3). We confirm the presence of stellar matter around the primary. Even after 
subtraction of ADS 9706B influence, light curve solutions show third light in the system. Models in FOTEL and Phoebe are presented. 
Conclusions. TW Dra is an astrophysically very interesting eclipsing binary. Future combination of interferometry, spectroscopy, and 
photometry is promising. It could definitely confirm the hypothesis of quadruple system for TW Dra and explain behaviour of this 
system in complex. 

Key words, stars: binaries: eclipsing - stars: fundamental parameters ~ stars: individual: TW Dra 



- ' — \ 

X 



(N 

> 

Qv^ 1. Introduction Although TW Dra was studied spectroscopically several 

f — \ times since 1919, secondary component was firstly detected in 

^ ' TW Draconis («=15^33-5P.l, <5=63°54'26"(2000.0)) is a well- ^P^^'^'^ by Smith ( 1949) and the radial velocity of both com- 

lO 1 J f» u J A1 1 » 1- • u- T-u • ponents were measured only by Popper ( 1989). The published 

: known and often observed Algol-type eclipsmg bmary. The vari- ^, ^ ,r^,-^^nT,~i 

(N ■ ability of the star was discovered by Cannon (see (Pickering °f *^ 'T^T''",'..'. ""^^ 

;[i9i^'. The variable star is an A-component of ^iiil^ EMg) to 0.47 Richards & Albnght Om)- 
. nary ADS 9706. The main light changes with amplitude in 

^,]B~ 2.3 mag in primary minimum ai-e caused by 11.5 hours TW Dra was also included in several surveys. 'Singh et alj 

> , eclipses of the hot main sequence A8V star by the cooler and ([T995I) studied X-ray radiation in five Algol-like stars using 

fainter KOIII giant component. ASCA and ROSAT satellites. He obtained first X-ray spek- 

trum of TW Dra and confirmed previous paper revealin g 

H First hghtcurvesolution TW Draconi s as a X-ray source (IWhite & Marshall [T983h . 

. ^ was made by ^aglowl (|I952D, who also mentioned large distor- |Umanaetal.l dMI) studied microwave radiation among Algols, 

tion of the secondary component. However, to obtained good r^^^ ^^^^ radio-fluxes of Algols ai-e generally comparable to 

solutions with low residuals he assumed very low values of ^^^^ ^^^^^^^^ ^^^^ CVn stars. During this study in the 

hmb darkening coefficients and he obtained also surprisingly -^^^^^^^ 1984-1989 went the radio-flux of TW Draconis down 

low reflection pai-ameters. During further and longer monitor- f^^^^^ summary although TW Draconis is known and 

ing of TW Dra behaviour Iwytd (11978D found large changes ^^^^^^^ ^^e satisfactory complex s olution offl iis unique system 

of orbital period and changes of the hght curve outside the j^^^ ^^^^ published up to now. IZeida et alJ JM) published 

minima, especially around secondary minimum caused by gas ^^^^-^^^ ^^^^^ ^^^-^^ ^^^^^^^^ j„ ^^^^^^ 
stream and a large hot spot in the target area of it. Observations 
made by Baglow (1952); Walter (1978) were reproce s sed b y 

e.g. Giuricin et al. (1980) or Al-Naimiv & Al-SikabI (11984 . The main goal of this study is to find the correct spectro- 

[Papousek et al. ( 1984) made several years to obtain light curve scopic mass ratio of the main components and to calculate new 

solution of the system. They mentioned larger scatter of the light curve and radial velocity curve solutions based on our new 

phased light curved but without any discussion in detail. measurements. 
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Table 1. Journal of photometric observations: Columns: 1. 
Observer(s); observatory. 2. Epoch. 3. Diameter/focal length of 
used telescopes in mm. 4-8. Numbers of measurements in pho- 
tometric bands. 



Observer(s); place 



Yeai"(s) Equipment(s) U B 



V R 



CCD 

Brat; Pec pod Snezkou 2005 
privat observatory, CZ 
Hroch; Brno, Masaryk 
university Observ.,CZ 
Chrastina, Szasz; obs. 2005 
Hlohovec, Slovaitia 
Svoboda; Brno, 2006-7 
privat observatory,CZ 
Zejda; N. Copernicus 2001-7 
Obs. and planet. Brno, CZ 



200/1850 
2006 620/2780 



179/1000 

15/116 
35/135 

400/1750 
40/360 
40/360 

200/1000 



154 - 
57 57 - 

300 300 299 301 
147 - 

688 - - 513 

225 6213 6089 6073 

900 3301 

1662 639 715 - 

7452 5017 5028 5021 



PEP 

Zejda, Janfk, Bozic; 2005 
Hvar observatory, Croatia 
Zejda. Janfk, Ogloza; 2006 
Mt. Suhora obs., Poland 



650/7280 136 136 136 - 
600/7500 - 2547 2552 2547 406 



2. Observations and reduction 

2.1. Photometry 

TW Dra was monitored during photometric campaign in 2001- 
2007. Almost 50000 usable measurements in UBVRI bands 
were collected (see journal of observations in Tab.lTJ. However, 
only Hvar measurements were transformed into the standard 
magnitudes. Suhora observations are only in differential mag- 
nitudes as well as all CCD measurements. 

The period of t he system is changing in at least two modes 
dZeida et al.L l2008l) . The present photometric data were all ob- 
tained in the epoch after the last great period change. Thus the 
period in the epoch of our observations is almost constant and 
light ephemeris in a linear approximation for this epoch could 
be used. Mikulasek's meth od based on Principal Component 
Analysis (hereafter PCA) (Mikulasek, 2007) was used to com- 
bine all our photometric measurements on different magnitude 
scale made in different time intervals to determine following 
light ephemeris: 

Pri.Min. = HJD 2453558.91888(32) + 2^806832 11 (28)-£. (1) 

Despite of our CCD observations lower accuracy, small os- 
cillations were detected on light curve, in B and V filters, respec- 
tively. Pulsati on period . 0519(^ 3) d and semi-amplitude Am - 
10(3) mmag (IZeida et al.L 120061) found in our B measurements 
and 0.0501(1 3)d with Am - 9(2) mmag in B measurements 

^rpn Z7~m Ai r\o yiK :„ *. ...:iU ry 



of [P apousek et al.l d 19841) are in agreement with Kusakin et aT 
(l20bl.) (0.0556 d. Am = 2.1 mmag in B) and Kim et al., (,20030 
(0.053 d. Am - 5 mmag). Change in pulsations amplitude is 
probably due to super position of diff erent pulsation modes as 
described in RZ Cas ( Mkrtichian et al.. 2008) . Large scatter in 
maximum of brightness in eclipsing phased Ught curve is also 
caused by oscillations and their superpositions from different or- 
bital cycles. 

Regardless of confirmed oscillations, deformations of light 
curve, which persist at least during several orbital cycles, can be 
easily seen as well. Our measurements show a small decrease 
(0.04 mag in B) after the egress from primary eclipse (phase 
0.10-0 .19) previously mentioned by Baglow ( 1952) and Walter 
( 119781) and a bump of 0.055 mag (V), 0X)30 mag (R), 0.02 mag 
(/) before the secondary eclipse (phase 0.40-0.43) (see Fig.[TJ. 



The most obvious light curve phenomenon is a dependence 
of the secondary minimum depth on wavelength. While in B 
band the secondary minimum is scarcely visible (Am=0.063 
mag), in / filter is easily detectable (Am=0.196 mag) (see Fig. 

m. 
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Fig. 1. Normal light curves around secondary minimum. 



2.2. Spectroscopy 

The spectroscopic observations (46 spectra) were made in JD 
interval 24531 1 1-2454218 using Coude spectrograph of the 2.0- 
m Ondi'ejov reflector, equipped with a SITe-005 800x2000 CCD; 
the spectra cover the Ha region from 628.0 to 672.0 nm, with a 
linear dispersion of 17 nmmm"', two-pixel resolution of 12700, 
and S /N ranging up to 478. Spectra with an average signal-to- 
noise (S/N) ratio lower than 95 in continuum near Ha line were 
excluded. Finally 37 spectra were analysed (see journal of our 
spectroscopic observations in Tab.|2]i. 

In all cases, the wavelength calibration was based on ThAr 
comparison spectra taken before and after the stellar exposure. 
Mean flatfield images from the same night as the stellar spectrum 
were always applied. 

The initial reductions to 1-D frames were carried out in IRAF 
by MS. The final reduction, rectification were carried out in 
SPEFO (iHorn et al.Lfl996LlSkodaL 119961) . 



3. Radial velocities 

We used three methods to obtain RV curves of components of 
this binary system. Unfortunately the lines of the secondary 
component being poorly resolved. The secondary component 
was detectable only on some well exposed spectrograms. In 
SPEFO the RV's were measured interactively, comparing the di- 
rect and flipped images of the line profiles. The zero point of RV 
scale was correcte d through the use of reliable telluric lines - see 
iHorn et all (Il996h for details. 

To obtain RVs we also used the cod e based on cross- 
corelation method (see IZverko et al.L l2007h . The needed syn- 
thetic S2e£tta__were_cal£ulated_accor to the parameters 
from lAl-Naimiv & Al-SikabI (11984 by the code SYNTHE (see 
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Fig. 2. Spectrum of TW Dra system in the phase 0.509 in com- 
parison with synthetic spectra of primary component (upper pan- 
els). Observed TW Dra spectrum in phase 0.001 (containing 
secondary component only) as compared with synthetic spec- 
trum of secondary component. Synthetic spectra are based on 
Kuru cz's model of atmospheres with following parameters taken 
from lAl-Naimiv & Al-SikabI (fT98l : T^g^i = 8300 K, log^i = 
3.92, visin; = 50 km s"', Teft^a = 4800 K, loggz = 3.24, 
V2 sin i = 72 km s"', solar abundance both. It is apparent at ob- 
served secondary component spectrum that absorption line Ho- 
is here dimmer than in synthetic spectrum, probably due to an 
emission in the system. 




643.0 e44.D e45.a etsD G47.D Bisa e4^.a e50.a 

wwrBlongEh [nm| 

Fig. 3. KOREL disentangUng - region 642.6-65 1 .0 nm. 



eral sources using Fourier disentangling with intrinsic line- 
profile variations. Twelve short regions in wavelength interval 
626-677 nm and later on 5 regions from them joint in one mul- 
tiple solution were used (see Tab.[3]l. 

We first measured the S /N for all spectra in the region 
661.5-662.0 nm and weighted them with weights proportional 
to (S /N)^ and normalized to one for the mean value. These 
weights we re used in the co ntrol file of the auxiliary program 
PREKOR dHadraval l2004bl) that prepares the input data for 
KOREL. Firstly we derived the line strengths of the telluric 
spectra for the wavelength range 646.5 - 651.0 nm, which con- 
tains a number of strong water vapour lines and then kept these 
strengths fixed in all subsequent solutions with KOREL. We de- 
rived the KOREL solutions separately for the regions shown in 
the Tab. [3] See an example of done disentagled region in Fig. [3] 

In the process of trial solutions was found that in wave- 
length interval 650-662 nm the Ha line profile is deformed by 
presence of emission (see Fig. |4]i, probably caused by the ac- 
cretion disc, which i s in ag reements with preliminary result in 
iRichards & Albright! (Il999l) . However, this region was not in- 
cluded into following solutions. 




ISbordone etaTl l2004lPl and the cod e SYNSPEC, which is able 
to calculate spectra from TLUSTY jLanz & Hubenvl l2007h and 
ATLAS model atmospheres. 

Results of the two first methods are shown in Table |2] 
Finally we used also the spectral disentangling technique (see 
relevant section) to obtain radial velocities. 

3.1. Spectral disentangling 

As a third approach for determin ation of radial velocit y 
curves the code KOREL was used (iHadravai 1 19971 l2004bh . 
This code enables decomposition of light joint from sev- 



* For the most recent update of the program and its manual 
see Sbordone & Bonifacio 2005, http:// wwwuser.oat.ts.astro.1t/atmos/l 

atlas ji'ookhook/ Atlas_Cookhook.html 




pnmary 

secondary 
accretion disc 
teluric lines 



Fig. 4. KOREL disentangling - Ha region. The emission com- 
ponent of spectra was disentangled as co-rotating object close to 
primary star 
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Table 2. Radial velocities in km/s measured in SPEFO (column 6-7), obtained from CCF (column 8-9) 



No. JD,„i Phase S/N Weight SVj KVi RVj RV2 phase 



1 


53111.5121 


0.6008 


162.4 


0.365 


35.4(1.0) 


-90(4) 


34.13(3) 




0.6008 


2 


53145.3644 


0.6615 


227.1 


0.713 


53.8(1.8) 


-106(3) 


56.80(3) 


-129.33(49) 


0.6615 


3 


53153.3578 


0.5093 


489.3 


3.311 


4.8(1.2) 


2(3) 


5.67(3) 


3.06(30) 


0.5093 


4 


53254.3689 


0.4969 


311.6 


1.343 


1.1(0.9) 


-4(4) 


3.10(3) 


7.73(32) 


0.4969 


5 


53454.4151 


0.7681 


389.0 


2.093 


66.0(1.3) 


-103(5) 


63.05(3) 


-134.91(49) 


0.7681 


6 


53459.4299 


0.5547 


98.8 


0.135 


16.9(1.3) 


-61(10) 


19.65(4) 




0.5547 


7 


53460.3532 


0.8837 


195.6 


0.529 


42.3(1.7) 


-98(4) 


41.46(3) 




0.8837 


8 


53461.5080 


0.2951 


219.9 


0.669 


-58.7(1.5) 


131(8) 


-60.71(3) 


133.83(51) 


0.2951 


9 


53463.4444 


0.9850 


106.8 


0.158 


-1.5(2.7) 


-14(4) 


9.94(19) 


-7.24(07) 


0.9850 


10 


53463.4912 


0.0016 


359.4 


1.786 


2.2(1.1) 


2(2) 


4.50(29) 


2.29(07) 


0.0016 


11 


53464.4606 


0.3470 


304.5 


1.282 


-50.7(1.4) 


108(13) 


-52.39(3) 


99.43(50) 


0.3470 


12 


53465.5142 


0.7224 


224.6 


0.698 


61.8(1.5) 


-110(6) 


65.20(3) 


-134.91(47) 


0.7224 


13 


53504.4166 


0.5823 


262.5 


0.953 


28.1(1.0) 


-101(5) 


27.75(3) 




0.5823 


14 


53510.4323 


0.7255 


305.8 


1.293 


59.3(1.5) 


-108(5) 


62.24(3) 


-134.91(53) 


0.7255 


15 


53511.4030 


0.0713 


274.5 


1.042 


-23.1(2.6) 


38(7) 


-23.02(4) 


-4.41(28) 


0.0713 


16 


53516.5391 


0.9012 


183.3 


0.465 


37.9(1.8) 


-97(4) 


43.22(4) 




0.9012 


17 


53538.3802 


0.6826 


243.9 


0.823 


57.2(1.8) 


-107(3) 


59.18(4) 


-129.42(48) 


0.6826 


IS 


53542.4899 


0.1468 


270.9 


1.015 


-44.7(1.8) 


102(7) 


-49.02(4) 


120 


0.1468 


19 


53550.4130 


0.9696 


137.5 


0.261 


11.8(3.7) 


-9(4) 


30.08(6) 


2.56(19) 


0.9696 


20 


53550.4357 


0.9777 


114.9 


0.183 


4.1(2.7) 


-13(4) 


17.54(9) 


-4.63(13) 


0.9777 


21 


53619.3175 


0.5184 


386.4 


2.065 


7.7(1.3) 


-10(4) 


9.08(3) 




0.5184 


22 


53764.6405 


0.2932 


189.4 


0.496 


-61.3(2.0) 


129(5) 


-55.49(3) 


146.73(52) 


0.2932 


23 


53764.6771 


0.3062 


223.0 


0.688 


-54.2(1.4) 


130(8) 


-54.14(3) 


132.75(55) 


0.3062 


24 


53765.6086 


0.6381 


372.3 


1.917 


45.7(1.2) 


-102(4) 


49.51(3) 


-93.09(51) 


0.6381 


25 


53846.3238 


0.3948 


351.6 


1.710 


-35.8(0.9) 


109(4) 


-36.33(3) 




0.3948 


26 


53860.3570 


0.3944 


364.9 


1.841 


-34.6(1.0) 


103(7) 


-33.91(3) 




0.3944 


27 


53991.3876 


0.0772 


196.2 


0.532 


-27.3(1.3) 


61(3) 


-24.87(4) 


-7.60(29) 


0.0772 


28 


54000.4012 


0.2885 


330.3 


1.509 


-62.0(1.4) 


144(5) 


-64.17(4) 


140.44(51) 


0.2885 


29 


54017.3270 


0.3187 


175.8 


0.427 


-56.3(1.1) 


118(7) 


-59.18(4) 


125.60(52) 


0.3187 


30 


54026.3265 


0.5250 


211.0 


0.616 


8.5(1.2) 


-8(4) 


9.94(3) 


3.10(28) 


0.5250 


31 


54115.6199 


0.3378 


187.9 


0.488 


-45.4(1.7) 


94(11) 


-45.37(3) 


97.30(49) 


0.3378 


32 


54191.4998 


0.3718 


149.7 


0.310 


-37.8(2.7) 


78(11) 


-40.29(3) 




0.3718 


33 


54192.4449 


0.7086 


319.6 


1.413 


61.7(2.6) 


-94(5) 


64.13(3) 


-142.46(48) 


0.7086 


34 


54192.4913 


0.7251 


191.9 


0.509 


58.7(2.9) 


-100(8) 


64.76(3) 


-144.31(49) 


0.7251 


35 


54203.3567 


0.5961 


364.4 


1.836 


39.8(1.3) 


-101(3) 


38.76(3) 




0.5961 


36 


54217.3376 


0.5772 


214.2 


0.634 


31.6(1.1) 


-92(4) 


30.62(3) 




0.5772 


37 


54218.4308 


0.9666 


254.2 


0.894 


29.0(2.1) 


-36(6) 


27.34(5) 


-0.49(21) 


0.9666 



Notes: Phases were calculated using ephemeiisfTl 



KOREL was able to overcome problems with faint and 
blended lines of secondary star. Nevertheless, first solutions (in 
resolution 2048 bins) gave two mass ratio values with compa- 
rable accurac y; qi - 0.465 close to the result of SPEFO, CCF 
an d result of Poppen ( 1989) and qo - 0.408 close to solution 
in iLehmann et al .1 (I2009h . This ambiguity can be found in pre- 
vious ly published valu es of q, which differ f rom 0.28 (iPearce , 
19371) or 0.37 ( Plaskett . [19191 iHohnberd Il9l to 0.465 (iPoppei , 
1989h or 0.47 (iRichards & Albrightlll994 



Using region No. 6 (see Tab. |3]l we built a grid of solution 
in resolution 512 bins for input values of q varying from 0.370 
to 0.528 with a step 0.002 and semi-amplitude of primary star 
radial velocity K\ in interval 61.0 - 66.8 km/s with a step 0.2 
km/s. During each solution q and K\ parameters were fixed. We 
obtained 2400 solutions in total. Analogous calculations were 
repeated for resolutions 256, 512, 1024, 2048 (corresponding to 
the widths 15.3, 7.6, 3.8 and 1.9 km/s/bin) with smaller steps of 
input values: q in 0.380 - 0.500, step 0.001; Ki in 60.0 - 66.0, 
step 0.1. Thus we obtained 7381 solutions for each resolutions. 
The global minimum in the space of solutions was found around 
q X 0.40 and thus found the final solution for all chosen spectral 
intervals (see Fig. |5j. The resulting mass ratio is ^ = 0.405(3). 
The resulting values for the best solutions are summarised in 
Tab. |4] Note that disentangling cannot provide information on 
the systemic velocity. Its value can only be obtained by a direct 
measurement of disentangled line profiles. This way, we found 
that the systemic velocity of the system is about -2 km s"'. The 
final radial velocity curves could be compared with other ones 
in Fig.|6] 




Fig. 5. Grid of solutions for varying entrance values q and Ki in 
KOREL, wavelength interval 642.6-651.0 nm. The global min- 
imum for best solution is the lowest value and the darkest 
place in the diagram. 
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Table 3. Selected spectral regions used for disentangling. 



200 



No. Region 



Primary comp. 



Secondary comp. 





[nm] 


line 


A 

[nm] 


line 


A 

[nm] 


1* 


632.50-636.97 


Si 11 


634.7109 


Cal 


633.0850 










Fe 1 


636.2338 


2 


633.97-642.69 


Si 11 


634.7109 


Fe 1 


636.2338 






Si 11 


637.1371 


Fel 


638.0743 






Fel 


639.3601 


Fel 


641.4980 






Fel 


640.0001 










Fel 


640.8018 










Fel 


641.1649 






3 


634.90-642.76 


Si 11 


637.1371 


Fe 1 


636.2338 






Fel 


639.3601 


Fe 1 


638.0743 






Fel 


640.0001 


Fel 


641.4980 






Fel 


640.8018 










Fel 


641.1649 






4* 


636.85-641.35 


Si 11 


637.1371 


Fel 


638.0743 






Fel 


639.3601 










Fel 


640.0001 










Fel 


640.8018 










Fel 


641.1649 






5 


641.50-646.03 


Cal 


643.9075 


Cal 


644.9808 






Fell 


645.6383 






6* 


642.60-651.00 


Cal 


643.9075 


Cal 


644.9808 






Fell 


645.6383 


V 1 


647.1662 






Cal 


646.2567 


Fe 1 


649.4980 


7 


646.50-651.07 






V 1 


647.1662 










Fe 1 


649.4980 


8 


650.50-662.96 


Ha 


656.2817 


Til 


662.5022 






Fel 


659.2914 






9 


661.40-666.07 






Ti 1 


662.5022 










Fe 11 


663.3750 










Sil 


664.3629 


10* 


664.00-669.69 


Fel 


667.7987 


Sil 


664.3629 






Hel 


667.8151 


Fe 1 


666.3442 


11* 


665.50-670.20 


Fel 


667.7987 


Fe 1 


666.3442 






Hel 


667.8151 






12 


667.00-671.11 


Fel 


667.7987 










Hel 


667.8151 







Notes: The asterisks in the first column shows if the corresponding spectral region was 
included in the joint 5-region solution. 



Table 4. KOREL solutions for the system for selected spectral 
regions. All epochs are given in HJD-2400000. Period is fixed 
P = 2'?80683211. 



Reg 


'on T p^ij^^tr. 




K, 


q 




1 


53558.9323 


64.1 


160.7 


0.399 


0.1996 


2 


53558.9342 


62.7 


155.4 


0.403 


0.1865 


3 


53558.9347 


62.3 


158.5 


0.393 


0.1771 


4 


53558.9335 


63.1 


160.6 


0.393 


0.1987 


5 


53558.9374 


63.9 


151.8 


0.421 


0.2212 


6 


53558.9239 


62.6 


154.7 


0.405 


0.2191 


7 


53558.9207 


62.8 


158.5 


0.396 


0.2009 


9 


53558.9157 


63.5 


156.1 


0.407 


0.1663 


10 


53558.9285 


63.7 


154.2 


0.413 


0.1747 


11 


53558.9212 


63.3 


157.1 


0.403 


0.1723 


12 


53558.9258 


66.2 


157.6 


0.420 


0.1555 


5 re 


gions 53558.9219 


62.7 


153.7 


0.408 


0.2168 



3.2. Initial analysis of RVs 

We also compiled all RVs available in the astronomical litera- 
ture and used the HEC19 prograrrQ to derive the heliocentric 
Julian dates (HJDs hereafter) for them. In particular, we used 



^ This program which handles data of various forms, together with 
brief instructions on how to use it, is freely available on the anonymous 
ftp server at http://ftp.astro.troja.mff.cuni.cz:hec/HEC19,. 




0,6 0.8 
phase 

Fig. 6. Radial velocities determined by SP EFO, CCF an d 
KOREL in comparison with values taken from Popper ("1989|). 
Phase was calculated according to the found mathematical 
model of period development. For more detail see IZeida et alJ 
(l2008l) . 

Table 5. Journal of available spectral observations and RV mea- 
surements of TW Dra. 



Obs.* 


Time interval 


No. of 


Source 




(HJD-2400000) 


RVs. 




1 


22062.9-22157.7 


14 


Plaskett (1919) 


2 


32342.6-32377.7 


55 


Smith (1949) 


3 


40290.0-45779.0 


77 


Popper (1989) 


4 


49077.8-49103.8 


104 


Richards & Albright (1999) 


5 


49108.7-49139.6 


130 


Richards & Albriaht (1999) 


6 


53335.4-53920.4 


37 


this paper 



*) In column "Ob.s'.", individual instruments are identified by numbers: 1... Dominion 
Astrophysical Observatory 1.83-m reflector, pg plates; 2...2-m McDonald telescope, 
McDonald observatory, Cassegrain spectrograph, pg plates; 3...3-m Shane telescope at the 
Lick observatory, pg plates; 4... 1.5 m National Solar Observatory (NSO) McMath-Pierce 
Main Telescope at Kitt Peak, echelle CCD spg.; 5... 0.9 m Coude Feed Telescope at Kitt 
Peak National Observatory (KPNO), CCD spg.; 6... Ondi'ejov 2-m telescope, coude CCD 
spg- 



the early photographic radial velocities from the DAO (^PlasketJ, 
1919), McDonald observ atorv (.Smith. .1949) and Lick observa- 
tory obtained by Poppej (Il989h . A journal of available spectral 
observations is in Table|5j 

All tr ial solutions of the RV were derived with the FOTEL 
program (lHadravaL[T990ll2004ah : for the results see Tab. |6] 



4. B-component of visual binary ADS9706B 

TW Draconis is A-component of visual binary OS 299 = ADS 
9706. The B-component is the star HD 140512=TYC 4184 61 2 
9.987 mag (VT) only 3.3" far away from TW Dra. ADS9706B 
itself has been photometrically observed only rarely. The most 
accurate measurements were done by Hipparcos satellite, where 
A-component (TW Dra ) is Hp = 7.517 + 0.006 mag and 
ADS9706B Hp = 9.887 ± 0.047 mag. Both components were 
measured together during our photometric measurements. To 
subtract the influence of B-component we obtained also three 
spectroscopic observations of this star. Two Ha region spectra of 
B-component from Ondfejov observatory and one echelle spec- 
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Table 6. Journal of RVs curves solutions using FOTEL. 



Ref. 


P [day] 


To (JD-2400000) 


Ki [km/s] 


q 


1 


2.80654 


22139.605(14) 


65.1(16) 




2 


2.806 


32341.891(17) 


62.9(21) 




3 


2.80686656 


42258.483(9) 


63.0(6) 


0.465(7) 


4 


2.80686656 


42258.501(8) 


63.3(6) 




5 


2.80683211 


53558.916(12) 


60.9(10) 


0.474(18) 


6 


2.80683211 


53558.938(16) 


62.7(9) 


0.468(23) 


7 


2.80683211 


53558.938(8) 


62.6(12) 


0.441(13) 


8 


2.80683211 


53558.9219(15) 


62.7(6) 


0.405(3) 



Notes: In col umn ''Ref \ the RVs sourc es are identi fied by numbers: 1.. Plaskett 1 1919), 2... 
ISmitHiT94g|) . 3...[Ppppei i 1989), 4... Poppeil iT98g|) bv observations, 5... this paper, SPEFO 
method, 6... this paper, CCF method, 7... this paper, CCF method, second template, 8... this 
paper, KOREL. 



trum from Tautenburg (iMkrtichian & Lehmanni l2"007h were ob- 
tained. Using equivalent width of chosen spectral lines and com- 
parison to the synthetic spectra, the spectral type of ADS9706B 
as F5V-G0V was determined. Measured Hp magnitudes were 
transformed into V magnitudes in Johnson system, according to 
[Harmanec ( 1998). The result for GOV spectral type was closer to 
our observational data. Thus, this spectral type and consequent 
parameters were used to subtract ADS9706B influence from our 
TWDra photometry. 

Two-colours diagram was created using decomposed bright- 
nesses of both components and colour indices of stars in 
the close vicinity of TW Dra (taken from the photometric 



database dMermilliod et alil2007[)). Measured v alues we r e com 



pared to theoretical models from 



Golavl([T97l :IC rame 
Considering this diagram, small distance of TW Dra d - 122 + 
15 pc (Ferryman & ESA, 1997) and high galactic lat itude (45°), 
interst ellar extinction was omitted. According to IChamblissI 
(Il992h . visual B-component could be a physical component of 
the system (along with TW Dra ). Using spectral type (GOV) 
of ADS9706B, distance 118 pc was obtained. This value is in 
agreement with the distance of TW Dra itself and thus supports 
Chambliss's conclusion. Though first astrometric measurements 
of ADS9706 were done in 1843, there is no significant change 
of position angle or angular distance up to now. 



5. Combined light-curve and orbital solution 

For full- fledged modelUng of all available data we used the code 
FOTEL ( Hadrava, '2004a') and the PHOEBE program version 
0.29c (Prsa & Zwitter, 2005; Prsa, 2006) with the WD back-end 
(I Wilson & Devinnevlll971l) . To improve the S /N of photomet- 
ric data, we created 1000 normal points in BVRI colors; in U 
band we used all 136 original individual measurements. We ex- 
cluded all observations from obtained with higher scatter in case 
of either a technical problem or bad weather RVs for the final 
solutions were taken from the KOREL solution. We fitted the 
model to photometric and RV data simultaneously. The limb- 
dark ening coefficients fo r the square-root law were interpolated 
from lvan Hammd d 19931) tables. 



5.1. FOTEL 



Photometric measurements of TW Dra from i Baglowl (Il952h : 
IWalteH (119781) ; iPapousek et all ([1984) were resolved as well. It 
is necessary to consider third light in the system to improve to 



— — 



"A r 




Fig. 7. Light curves of TW Dra together with solution from 
PHOEBE in UBVRI bands (from upper left). Radial velocity 
curve with the solution from PHOEBE (bottom right). 



quality and accuracy of results. Presence of the third light is in 
agreement with previously published solutions. 

Our own measurements were split into two datasets - photo- 
electric photometry from Hvar observatory transformed into the 
standard UBV magnitude and all (mainly CCD) observations in 
differential UBVRI magnitudes. Hvar measurements are approx- 
imately the same amount as Baglow's ones, but unfortunately 
not covering secondary minimum. Thus, solution with the third 
light fixed was calculated in this case only. Initial light ephemeris 
is: P=2.80683211 d, Mq = 53558.91888. The period P, mass 
ratio q - 0.405 and the eccentricity e = were fixed. Initial val- 
ues of t h e othe r parameters {r\,r2,i,q,K\) were adopted from 
iBaglowl (Il952h . Limb d arkening coefficien ts were interpolated 
according to tables from va n Hammd (11993 ). 

The second, CCD datasets include normal light curves in 
BVRI bands, each of them with 1000 points. Hereafter only these 
normal CCD light curves were utilised, however initial param- 
eters were the same as described above. Resulting parameters 
are given in the first column of Tab. |7] Solutions based on both 
datasets are in very good agreement. 

5.2. PHOEBE 

According to experience with our datasets in FOTEL we have 
proceeded solutions with non zero third light in PHOEBE only. 
First, radial velocities were used to improve initial values of pa- 
rameters associated with spectroscopy with following results: 
semi-major axis a - 12.0124(3) Rq, mass ratio q - 0.4026(11), 
inclination / = 87.13(5)°, system velocity 7=-0. 713(23) km s"'. 
The a, q, y parameters were then fixed. 

Input parameters were the same for both PHOEBE and 
FOTEL analysis, except above mentioned a,g',!,7 parameters. 
In contrast to FOTEL, temperatures of components were fitted 
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Table 7. Table of TW Dra light curve solution based on photo- 
metric data obtained during the campaign 2001-2007 (all with 
non zero third light). Period was fixed P = 2.8068321 1 d. Fixed 
parameters are marked b y an asterisk. Full version of the table is 
available in lZeidal (1200 8i) . 





band 


FOTEL 


PHOEBE 






including 


mode 2 


mode 5 






3rd hght 


(detached) 


(algols) 


To 










(53558. H-) 




.91894(4) 


.92013(25) 


.91879(6) 


ri 




0.2166(5) 


0.2139 


0.2106 


^ I, pole 




0.2151 


0.2127 


0.2096(14) 


^ I, point 




0.2180 


0.2157 


0.2124(17) 


^ I, side 




0.2166 


0.2142 


0.2110(15) 


f I, hack 






0.2153 


0.2120(16) 


ri 




0.3144(3) 


0.3047 


0.3030 


^2,pole 




0.2918 


0.2839 


0.2830 


^ 2,poini 




0.3529 




0.4077 


f2,side 




0.3077 


0.2961 


0.2950 








0.3294 


0.3277 


in 




86.57(6) 


87.30(8) 


87.13(3) 


a[Ro] 




12.084 










2.14 


2.112 


2.112(53) 


M2 [Mo] 




0.87 


0.850 


0.850(21) 






0.405 


0.4026* 


0.4026* 


7 




-0.7116 


-0.713* 


-0.713* 


Qi 






5.095(11) 


5.166(15) 


O2 






2.679(1) 


2.680* 


log^i 






3.940 


3.960 


log 82 






3.240 


3.240 


TdK] 






8180* 


8180* 


T2[K] 






4437(18) 


4407(23) 


Mboidmag] 






1.260 


1.260 


Mboiiimag] 






3.120 


3.160 


no-cf 


U 


0.0249 


0.0185 


0.0213 




B 


0.0140 


0.0102 


0.0099 




V 


0.0149 


0.0125 


0.0124 




R 


0.0138 


0.0122 


0.0109 




I 


0.0117 


0.0102 


0.0107 


S(0-C)2 RV 




0.8284 


1.7189 


1.7181 



in PHOEBE. Due to presence of third light it was impossible 
to use reliably either colour constraining method or the depth 
of minima for at least temperature ratio determination. Thus, 
temperature of primary star was e stimated according to spec- 
ti-al type A5V Ti = 8180 K dCoxl [l999.) . Based on the range 
of possi ble spectral t ypes of secondary component K2III, K OllI, 
K2V ( PoDPe?. ' 1 9891: 1 Yoon et al.L [19941: ISamus elail 12009^. cor- 
responding temperature of secondary star is in the interval 4390- 
4830 K (iCoxL I1999h . Solutions for the mode 2 (detached sys- 
tems) and mode 5 (algols) were proceeded and resulting param- 
eters are listed in 2nd and 3rd columns of Tab.|7l showing similar 
results with only small differences. Thus the secondary compo- 
nent fills its Roche lobe or is only slightly below the Roche limit. 

6. Conclusions 

Using rich photometric material we found new linear ephemer is 
(in agreement with previous period study (IZeida et al.L l2008l) ). 



Although this study was not aimed to detailed study of de- 
tailed photometric behaviour, we confirmed pulsation of one 
component in o ur data as well as in previous set od data from 
iPapousek et al.l ((19841) . There are also small irregularities in light 
curves caused by spot activity on surfaces of components or in 
surrounding stellar matter 

Spectroscopic part concentrated on several tasks. Radial ve- 
locities were obtained from new spectra applying three methods 
(SPEFO, CCF, KOREL); radial velocity curves for both compo- 
nents were created and mass ratio of component was found with 
the value q = 0.405(3) (KOREL, FOTEL) and q = 0.4026(11) 
(Phoebe solution), respectively. Using a disentangling technique 
a presence of interstellar material in the system presumably in 
a form of accretion disc was con firmed. Detailed study using 
set of old yet unprocessed spectra [Richards & Albrighf (1 19991) 
as well as new huge collection of spectra from Tautenburg 
(iLehmann et al.L l2008l) could help to study the presence of cir- 
cumstellar matter in the system and its evolution in time. 

The first obtained spectra of visual component ADS9706B of 
TW Dra served for its spectral type determination (GOV). This 
result supports previously published hypothesis that ADS9706B 
is not only visual but even a physical member of the system 
(IChambUssl |1992|) . However, IZeida et al.1 (l2008 f) found on the 
base on light time effect appeared in O-C residuals another com- 
ponent which is too faint to exert influence upon the light of 
two basic component and above mentioned third component 
(ADS9706B). Thus TW Dra itself is a basic binary in a quadru- 
ple system ((1-i-1)-i-1)h-1. The spectral type GOV also allows to 
subtract the influence of ADS9706B from total light of visual 
pair ADS9706. Nevertheless new more precise spectroscopic 
and photometric measurements of B -component are highly de- 
sirable. 

Radial velocity and light curve solutions were done in 
FOTEL using entire (historical and current) photometric and 
spectroscopic observations. New photometric and spectroscopic 
observations were processed in PHOEBE as well. The param- 
eters from both FOTEL and PHOEBE are in very good agree- 
ment. TW Dra belongs to small numbers of semidetached eclips- 
ing binaries with well known basic parameters. The nowadays 
primary components of these systems (e.g. KO Aql, S Cnc, RZ 
Cas, TV Cas, U CrB, Al Dra, S Equ, TZ Eri, AF Gem, TT Hya, 
BP Mus, AT Peg, HU Tau, TX UMa and others) obtaining ma- 
terial from the secondary present themselves like main sequence 
stars, while the the same age secondaries are bigger and more ra- 
diative. Obtained results of TW Dra correspond to other algols 
with known accurate parameters (see Fig.|8]l. 
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